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Abstract 

Flammulina velutipes is one of the major edible muslnrooms in the world. Recently, abnormalities that have a negative 
impact on crop production have been reported in this mushroom. These symptoms include slow vegetative growth, a 
compact mycelial mat, and few or even no fruiting bodies. The morphologies and fruiting capabilities of monokaryons of 
wild-type and degenerate strains that arose through arthrospore formation were investigated through test crossing. Only 
one monokaryotic group of the degenerate strains and its hybrid strains showed abnormal phenotypes. Because the 
monokaryotic arthrospore has the same nucleus as the parent strain, these results indicated that only one aberrant nucleus 
of the two nuclei in the degenerate strain was responsible for the degeneracy. A sequence-characterized amplified region 
marker that is linked to the degenerate monokaryon was identified based on a polymorphic sequence that was generated 
using random primers. Comparative analyses revealed the presence of a degenerate-specific genomic region in a telomere, 
which arose via the transfer of a genomic fragment harboring a putative helicase gene. Our findings have narrowed down 
the potential molecular targets responsible for this phenotype for future studies and have provided a marker for the 
detection of degenerate strains. 
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Introduction 

The basidiomycete Flammulina velutipes is one of the most 
cultivated edible mushrooms. One of the reasons for the recent 
increase in mushroom production is consumer interest in their 
high content of macromolecules that possess antitumor, immuno- 
modulatory, and antiviral properties, such as polysaccharides and 
glycoproteins [1,2]. 

Filamentous fungi that are grown on a nutritionally rich 
medium frequendy exhibit instability manifested in morphological 
and physiological variations. The morphological variations include 
a lack of or reduction in sporulation, fluffy mycelial-type growths, 
and variations in hyphal pigmentation [3,4]. Previous studies 
implicated infections with double-stranded RNA viruses [5] and 
the instability of fungal nuclear and mitochondrial genomes in 
changes in fungal morphology and physiology. Such changes 
might be expected because fungal genomes are dynamic and 
capable of a rapid accumulation of genome rearrangements, 
particularly amplifications and deletions [6]. Likewise, edible 
mushrooms that were preserved by serial passaging in culture on 
nutritionally rich media exhibited abnormal mycelial growth and 



poor yields [7] . Abnormalities in Agaricus hisporus led to a lower 
yield and inferior product quality, resulting in economic losses 
[8,9,10]. These abnormalities may be due to changes in the 
ribosomal-DNA copy number, the loss of heterozygosity at specific 
loci, de-heterokaryotization, chromosomal loss, or chromosomal- 
length polymorphisms [10]. In the case off. velutipes, abnormal- 
ities such as malformed fruiting bodies and the complete loss of 
fruiting-body development have been observed [11]. Despite the 
relative frequency of these occurrences, the underlying cause is not 
entirely clear. Notwithstanding its economic importance, many 
basic questions involving the biology of this fungus remain 
unanswered; hence, a method for detecting mutant strains of F. 
velutipes would be very useful. F. velutipes, unlike many 
basidiomycetes, produces abundant arthrospores in both mono- 
karyotic and dikaryotic myceha [12]. Moreover, most of the 
arthrospores contain only one haploid nucleus, which is generated 
without resorting to fruiting-body formation and chromosomal 
rearrangements [13]. Thus, the normality of each nucleus can be 
easily determined by test crossing. Additionally, molecular 
selection markers such as sequence-characterized amplified 
regions (SCARs) are required to assist mushroom farmers in 
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cultivating wild-type strains rather than degenerate clones. Herein, 
we describe the characterization and identification of a degenerate 
nucleus and report a specific marker for differentiating wild-type 
and degenerate nuclei. 

Materials and Methods 

Strains and pure cultures 

The wild-type (Fvl-5) strain and 3 degenerate strains (Fvl-5'", 
Fvl-5''^ and Fvl-5''^) of F. velutipes were obtained from various 
mushroom farms (Table 1). The mycelia were grown on 
mushroom complete medium (MCM) at 25°C. Because F. 
velutipes mycelia are typically transferred as a plug rather than a 
single cell, the collected degenerate strains most likely contained 
some normal cells. To enhance the purity of the degenerate 
strains, the mycelia were serially cultured on modified BTB- 
sawdust agar containing 0.45% peptone, 0.75% yeast extract, 
0.5% oak sawdust, 0.025% bromothymol blue (BTB, Sigma- 
Aldrich, St. Louis, MO, USA), and 2% agar, growth on which has 
been reported as an indicator of normal function [11]. At each 
round of isolation, the mutants that did not decolorize the BTB 
agar and exhibited abnormal morphologies, such as a reduced 
growth rate and compact mycelial colonies, were transferred to 
fresh medium via cells at the end of the mycelia. The isolated 
degenerate and wild-type strains were grown on MCM agar for 
subsequent analyses, and their mycelial growth rates along a 
perpendicular hue were determined at 3, 5, and 7 d (dikaryons) or 
at 8, 10, and 15 d (monokaryons) after inoculation. The 
experiments were performed in triplicate. 

Monokaryon isolation, test crossing, and fructification 

The purified wild-type and degenerate strains were grown on 
MCM agar for 7 d at 25°C. Arthrospores were obtained from the 
solid medium by adding 1 to 2 mL of sterile water to the cultures 
and gently prodding them for a few seconds. The suspension was 
serially diluted, plated onto MCM agar after filtering through glass 
wool, and grown at 25°C. The germinated mycelia were isolated, 
and the clamp connections were observed using phase-contrast 
microscopy at 40 x magnification (Olympus 1X71; Olympus, 
Tokyo, Japan). Monokaryons with no clamp connection were 



subjected to reciprocal crossing to identify their mating type 
(Table 1). To determine the abnormal monokaryon group in the 
two compatible groups that were derived from the degenerate 
strain, all possible test-crossing combinations were performed 
between monokaryons derived from the wild-type and degenerate 
strains. The hybrids were subjected to growth-rate analysis and the 
fruiting test. A substrate containing pine sawdust (23%), corncob 
(29%), rice bran (18%), beet pulp (4%), wheat bran (14%), 
cottonseed hull (4%), shell powder (4%), and soybean powder 
(4%), with a 65% water content, was autoclaved at 121°C for 
1.5 h, cooled to 20°C, inoculated with four plugs (1x1 cm) of 
cultured hybrids, and finally placed in an incubation room 
maintained at 19°C to 20°C for spawn running. On the 35th day 
of the spawn running period, the outer area of the substrate was 
removed by scraping to induce the formation of primordia. Then, 
the cultures were placed in a cultivation room that was maintained 
at 15°C and 90% relative humidity to induce fruiting. When the 
fruiting bodies reached a length of approximately 2 cm, the 
temperature was lowered to 5°C for 5 to 7 d to normahze the size 
of the fruiting bodies. Subsequently, the cultures were maintained 
at 10°C with 75 to 80% relative humidity (RH) to allow the 
fruiting bodies to grow. 

Extraction of genomic DNA 

Genomic DNA was extracted from dikaryons and monokaryons 
using the Exogene Plant SV kit (GeneAll Biotechnology, Korea) 
according to the manufacturer's protocol. The yields of DNA were 
quantified using a NanoDropND-1000 UV spectrophotometer 
(NanoDrop Technologies, Wilmington, DE, USA). 

Identification of an abnormality-associated marker 

A draft genomic sequence of the monokaryotic strain (Mono3 
derived from meiotic spore of Fvl-5) of F. velutipes that was 
obtained in a previous study (http://l 12.220.192.2/fve/) and the 
complete genomic sequence of KACC42780 [14] were used as 
sources for SSR marker development and comparative analysis. 
The SSR candidates were identified using SSR Locator I [15]. 
Primers were designed accordingly, using the Primer3 program 
[16]. Random-amplified polymorphic DNA (RAPD) PGR was 
performed in 20-|J.L mixtures containing 30 ng of genomic DNA, 



Table 1. List of the F. velutipes strains used in this study. 



strain 


Description 


Stock No* 


Fvl-5 


Wild-type, cultivar 


KCCM 90233 




Degenerate dikaryon 


KCCM 90234 




Degenerate dikaryon 




Fvl-S"^ 


Degenerate dikaryon 




Wl xW4 


Hybrid between monokaryons from a wild-type strain 


KCCM 90239 


Wl xD4 


Hybrid between monokaryons from wild-type and degenerate strains 


KCCM 90242 


Dl xW4 


Hybrid between monokaryons from wild-type and degenerate strains 


KCCM 90240 


Dl xD4 


Hybrid between monokaryons from a degenerate strain 


KCCM 90241 


Wl 


Monokaryon isolated from Fvl-5, AxBy 


KCCM 90235 


W4 


Monokaryon isolated from Fvl-5, AyBx 


KCCM 90236 


Dl 


Monokaryon isolated from Fvl-5''\ AxBy 


KCCM 90237 


D4 


Monokaryon isolated from Fv^-S'^^, AyBx 


KCCM 90238 



*: KCCM: Korean Culture Collection of Microorganisms. 
doi:l 0.1 371/journal.pone.Ol 07207.t001 
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Figure 1. Mycelial growth of wild-type and degenerate 
dikaryons of F. velutipes. (A) IVIorphology of the mycelia of the 
wild-type (left) and degenerate strains (right) that were grown on IVICIVI 
agar for 7 d. (B) Growth rates of the wild-type and degenerate strains 
were determined based on the diameter of the mycelia grown at 25 C. 
The error bars represent the mean standard deviations of triplicate 
samples. 

doi:1 0.1 371 /journal.pone.01 07207.g001 



1 X e-Taq buffer, 0.2 mM dNTPs, 0.2 |Xmol of Taq polymerase 
(Solgent, Korea), and 1 |a.M random primers of the OPB, OPC, 
OPD, OPE, OPF, OPG and OPI series (Operon Tech., 
California, USA). PCRs were performed using a Gene Atlas 
system (ASTEC, Japan) and the following protocol: an initial 
denaturation step of 5 min at 95°C; followed by 40 cycles of 1 min 
of denaturation at 95°C, 1 min of annealing at 40°C, and a 2-min 
extension at 72°C; with a final extension step of 5 min at 72°C. 
The SSR PGR cycles consisted of an initial denaturation step of 
3 min at 95°G; foUowed by 35 cycles of 20 s at 95°C, 40 s at 52°G, 
and 30 s at 72°C; with a fmal step of 5 min at 72°G. The PGR 
products were resolved using 1% (RAPD) and 3% (w/v: SSR) 
agarose gels (Life Technologies, USA) in TAE buffer (400 mM 
Tris, 200 mM sodium acetate, and 20 mM EDTA, pH 8.3) 
containing RedSafe (Intron, Korea). The specific degenerate DNA 
bands were excised from the gel, purified using an Expin PGR SV 
kit (GeneAU Biotechnology, Korea) as described in the manufac- 
turer's protocol, ligated into a vector using a Dr. TA TOPO 
cloning kit (Doctor Protein, Korea) and were then used for 
bacterial transformation. The recombinant plasmids were isolated 
using a Plasmid DNA purification Hybrid-Q kit (GeneAU 
Biotechnology, Korea) and were sequenced (Macrogen Gorp, 
Korea). To increase the reproducibility and reliability of the results 
obtained using the random primer set, we constructed a SGAR 
marker (22-mer) based on the sequence of the polymorphic RAPD 
band (Table SI). 



Comparative analysis 

The acquired sequences specific to the degenerate monokaryons 
were aligned with the genomic sequences of the monokaryotic 
strain Mono3 using the DNAMAN program (Lynnon Gorp, 
Ganada) to determine its flanking sequence. Because the 
degenerate-specific region of the degenerate (D4) strain was 
determined according to scaffold 59 of the Mono3 genomic 
sequence, PGR was performed using primer sets spaced every 1 kb 
(for the sequence from 1 to 10 kb), 10 kb (for the sequence from 10 
to 60 kb) or 60 kb (for the sequence from 60 to 121 kb) based on 
the scaffold 59 sequence for synteny analysis of the wild-type (W4) 
and D4 strains. To determine the genomic region corresponding 
to the degenerate-specific sequence and the boundary regions for 
the W4 and D4 genomic sequences, the primer sets for 9,244 to 
10,004 and 9,500 to 10,500 and other regions were used to 
amplify DNA fragments, which were sequenced (Table SI). PGR 
using the same mixture described above was performed in 20 |iL 
volumes using the following protocol: initial denaturation for 
5 min at 95°G, followed by 35 cycles of 1 min denaturation at 
94°G, 1 min annealing at 60°G, and a 90-s extension at 72°G. The 
PGR products were evaluated using electrophoresis on 1% (w/v) 
agarose gels. The PGR products were hgated into a vector using a 
Dr. TA TOPO cloning kit and sequenced in the forward and 
reverse directions using Ml 3 primers. The sequences were 
assembled using the DNAMAN program (Lynnon Gorp). Three 
sequences from the W4, D4 and KAGG42780 strains were aligned 
using the DNAMAN program. The sequence specific to the D4 
strain was analyzed to predict the presence of genes and the 
functions and domains of the gene products using the FGENESH 
program (http://www.softberry.com/), the UniProt program 
(http:/ /www. uniprot.org/) and the Protein BLAST program from 
the NGBL Tandem repeats in the genomic DNA were identified 
using the Tandem Repeats Finder program (http://tandem.bu. 
edu/trf/trfhtml). 

Results and Discussion 

Growth characteristics and morphological traits of the 
degenerate mycelium 

The degenerate strains exhibited abnormal phenotypes, includ- 
ing reduced growth rates and the formation of compact mycelial 
mats compared to the wild- type strain (Fig. lA); however, the 
extent of these abnormalities was strain-dependent (data not 
shown). The wild-type F. velutipes strain elicited a color change in 
BTB from blue to yellow, whereas the degenerate strains did not 
[11]. The results of the decolorization of the dye by these fungi 
suggest that laccase and peroxidase were responsible for the 
degradation activity via oxidative reactions [17]. The phenolic 
structure of BTB, which includes a methyl group, is similar to that 
of lignin and is targeted by ligninolytic enzymes just as lignin is 
degraded [18]. During the initial transfers (i.e., within the first and 
second rounds of transfer), the non-decolorized regions were 
scattered at the edge of the growing mycelial colony in the 
degenerate strains, whereas the wild-type strain completely 
decolorized the BTB. The decolorization activity was enhanced 
when the MGM + BTB was supplemented with sawdust (data not 
shown). After serial transfers of the non-decolorizing regions to 
medium supplemented with BTB and sawdust, the non-decolor- 
izing strains Fvl-5'", Fvl-5'*^, and Fvl-S"*^, which exhibited 
consistent phenotypes throughout the mycelial mats, were 
obtained. The growth rate of the degenerate mycelia (Fvl-5'^') 
was 1.87 cm/3 d, whereas that of the wild-type mycelia was 
3.36 cm/3 d (Fig. IB). Poor quality and yield have also been 
reported to be accompanied by the abnormal microscopic 
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Figure 2. Growth of monokaryotic mycelia and test crosses in F. velutipes. (A) Crossing table of the monokaryons derived from the wild-type 
or degenerate strains. (B) Morphology of the mycelia of the degenerate monokaryons (right) and their wild-type (left) counterparts. Wl and Dl (upper 
panel) and W4 and D4 (lower panel) have compatible relationships. (C) IVIycelial morphology of hybrids from crosses between monokaryons of wild- 
type and degenerate strains. The growth rates of the monokaryon (D) and hybrid test crosses (E) were determined based on the diameter of mycelia 
grown at 25 "C. 

doi:1 0.1 371 /journal.pone.01 07207.g002 



morphology of the myceha oiA. hisporus, including the branching 
pattern and thickened hyphae [8,9]. However, no difiFerence 
between the wild-type and the mutant strains was observed in our 
study, consistent with the results of the study of Magae et al. (2005) 
on mycelial morphology at the microscopic level (data not shown). 

To identify the abnormal-nucleus group, two compatible 
monokaryons of the wild-type strain and the degenerate strains 
were prepared through arthrospore formation followed by mating. 
Three isolates of each mating group of the wild-type and 
degenerate strains were selected for further study. The mating 
types were divided into two compatible categories, with the wild- 
type group "a" (Wl, W2, and W3) and the degenerate group "a^'" 
(Dl, D2, and D3) on one side and the wild-type group "b" (W4, 
W5, and W6) and the degenerate group "b^'" (D4, D5 and D6) on 
the other (Table 1, Fig. 2A). The growth rate of the mycelia in 
group b*^ was less than half that of its wild-type counterpart, but 
the mycelial mat of the former was more compact (Fig. 2B). In 
contrast, the morphology (Fig. 2B) and growth rate (Fig. 2D) of 
the other monokaryons (a'') were not distinguishable from those of 
the wild-type group (a). To elucidate the effect of monokaryotic 
abnormalities on the dikaryon, hybrids were obtained from inter- 
and intra-compatible combinations of the monokaryons of Fv 1 -5 
and Fvl-S"". Clamp connections were observed despite the 
presence of chromosomal abnormalities (data not shown), 
indicating that the formation of clamp connections is not 
influenced by the degeneration of these strains. The hybrids 



axb'' (Wl xD4, W2xD5, W3xD6) and a'^xb'* (Dl xD4, D2xD5, 
and D3xD6) exhibited phenotypes similar to those of Fvl-S"^', 
including slow propagation of the mycelia and the formation of a 
mycehal mat on the MCM agar plates (Figs. 2C, 2E). The 
phenotypes of the hybrids (axb and a''xb; Wl xW4, W2xW5, 
W3xW6, Dl xW4, D2xW5, and D3xW6) were not significandy 
diflFerent from those of the wild-type hybrids. Previous reports 
[19,20] have identified mitochondrial disorders as causative agents 
of mycelial abnormalities, but our results excluded that possibility 
because the group a"^ and b*^ monokaryons shared an identical 
mitochondrial pool, yet only one group showed a distincfly 
degenerate phenotype. 

Fruiting-body formation by the hybrids 

To compare the fruiting-body formation of the wild-type and 
degenerate monokaryons, hybrid mycelia were cultivated in 
bottles on sawdust medium at a commercial scale. The Dl xW4 
hybrid showed spawn running and developed characteristics 
similar to those of the wild-type strain (Fig. 3A). In contrast, for 
the Wl xD4 hybrid, the spawn did not completely run on the 
sawdust medium at 35 d nor were primordia or fruiting bodies 
observed (Fig. 3B). Other sibling combinations of Wl xD4 showed 
consistent results (data not shown). In contrast, the hybrids (the 
axb group and the a'^xb group) showed normal spawn running, 
primordial emergence, pinheading, and fruiting-body develop- 
ment that were undistinguishable from those of the wild-type 
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Figure 3. Mycelial running and formation of fruiting bodies. (A) Mycelial running of the wild-type and degenerate dikaryons and the hybrids 
of crosses between monokaryons of the wild-type (Fv1-5) and degenerate (Fv1-5^) strains grown on sawdust medium in bottle cultures for 35 d at 
20°C after inoculation. (B) Morphologies of the fruiting bodies produced by the strains represented in the left panel. The experiments were performed 
using triplicate samples. 
doi:1 0.1 371 /journal.pone.01 07207.g003 



strain (Fig. 3A, B). Considering the abovementioned observations, 
the b"^' group appeared to be responsible for the serious defect in 
fruiting-body differentiation. 
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Figure 4. PCR-amplification products of the RAPD and SCAR 
primers. (A) RAPD electrophoretic profile of monokaryons of the wild- 
type and degenerate strains obtained using random primer OPD-05. 
The arrows indicate the polymorphic bands of the wild-type (W4) and 
degenerate monokaryons (D4). (B) Specificity of the FEDE primer set 
(SCAR marker). Wild: Fv1-5 (wild-type), De: Fvl-5'' (degenerate). W1 to 
W6 are monokaryons from Fv1-5, and D1 to D6 are monokaryons from 
Fvl-5''. W1-W3: group a {AxBy), W4-W6: group b (AyBx), D1-D3: group a"" 
[AxBy), and D4-D6: group b'' {AyBx). 
doi:1 0.1 371/journal.pone.01 07207.g004 



Identification of an abnormality-associated marker 

Polymorphisms were frequently detected in the RAPD PGR 
products of compatible monokaryons (groups a'' and b''), whereas 
polymorphisms were very rare in the RAPD PGR products of 
groups b and b''. Of the 438 primers tested, only the OPD-05 
primer detected a polymorphism at approximately 2.1 kb 
(Fig. 4A). The SGAR marker was designed based on a specific 
amplified sequence obtained from a single fragment (Fig. 4B), 
which appears to be present in the b"^, a, and a"* strains but not in 
the b strain. However, this marker could not discriminate between 
the dikaryotic wild-type and dikaryotic degenerate strains. The 
SGAR-marker primer set consistently amplified the degenerate- 
specific band from other monokaryons (the same mating group as 
b) of the degenerate strains (Fvl-S''^ and Fvl-S'*'^) (data not shown). 
These results suggested that the transfer of a genomic region that 
included the SCAR primer binding sites from an a strain nucleus 
to a b strain nucleus. Although this screening procedure was not 
applied to other F. velutipes cultivars in the present study, FvI-5 is 
a widespread cultivar in Asia, and the degenerate strains collected 
from farms exhibited similar abnormalities, including compact 
myceha and slow growth. Thus, our procedure will be helpful to 
commercial farmers for distinguishing degenerate mycelia from 
the complex mixture of normal and degenerate structures that are 
frequently found on mushroom farms. 

Genomic localization of the region responsible for the 
abnormalities 

The alignment of the sequence acquired from D4 with the 
Mono3 genomic sequence [21] showed more than 99% sequence 
similarity with bases 2,859 to 5,041 of scaffold 59 (data not shown). 
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Figure 5. Synteny analysis of tKie genomic regions of D4 and 
W1 specific to tiie degenerate strain. (A) DNA fragment of the W4 
and D4 genomic region corresponding to 1 to 1 0 l<b of scaffold 59 from 
the F. velutipes monokaryotic strain (IVIonoS) obtained by PCR using the 
SCAR primer set (FEDE-F and -R). Wild: Fv1-5 (wild-type), De: Fvl-S"^ 
(degenerate). W1 to W6 are monokaryons from Fv1-5, and D1 to D6 are 
monokaryons from Fvl-5''. W1-W3: group a (AxBy), W4-W6: group b 
[AyBx], D1-D3: group a"* (AxBy), D4-D6: group b"" (AyBx). (B) DNA 
fragment of the W4 and D4 genomic region corresponding to 60 to 121 
kb of scaffold 59 obtained by PCR using the SCAR primer set. 
doi:1 0.1 371 /journal.pone.01 07207.g005 

Based on the synteny analysis of the wild-type and degenerate 
genomic regions corresponding to this scaffold, the region from 
approximately 1 to 10 kb showed that specific bands of the 
expected size could be detected in D4 but not in its counterpart 
W4 (Fig. 5A). In contrast, similar band patterns in the wild-type 
and degenerate strains were detected in the regions beyond 10 kb 
(Fig. 5B). To determine the boundary region, PCR was performed 
using the primers that hybridized at 9,244 bp in the genomic 
sequence that only yielded a band for D4. However, primers that 
hybridized at 9,500 bp yielded amplicons for D4 and W4, 
indicating that the boundary region of the wild-type and 
degenerate genomic region lie between 9,244 and 9,500 bp 
(Fig. 5 A, B). The alignment of the sequences from D4, W4 and 
KACC42780 showed that the region corresponding to the 
degenerated region (1 to 9,244 bp of scaffold 59) was detected 
only in the D4 genome, indicating the involvement of this region 
in the degenerate phenotypes (Fig. SI, S2). Comparative analysis 
of the genomic region between 9.5 to 9.7 kb of D4 with that of the 
KACC42780 genome revealed that the location of the degenerate- 
specific region is near the terminus of chromosome 8 (—116 kb 
from the end) (Fig. S2). At the 3' terminus of the D4 genome, 
multiple telomeric repeating units (TTAGGG) that were found in 
the telomeric region of P. ostreatus (Perez. 2009) were detected 
(Fig. SI), indicating that the degenerate-specific sequence is likely 
to be a telomeric region. The similarity between this region in 
scaffold 59 (Mono3) and its corresponding region in Ch8 
(KACC42780) ranged from 30 to 79% because Mono3 is a highly 
inbred cultivar, whereas KACC42780 was isolated as a wild-type 
strain. 

Based on the FGENESH, UniProt and protein blast analyses, 
the region corresponding to the degenerate-specific genomic 



region of the D4 strain encodes a protein of 2,432 amino acids 
and a molecular weight of 268 kDa, which showed similarity to 
the putative ATP-dependent helicase C17A2.12 of Thanatephorus 
cucumeris and contains HELICc and HepA domains (Fig. SI). 
These helicases are reportedly involved in various functions, such 
as replication, transcription, translation, recombination, DNA 
repair, and ribosomal biogenesis [22]. Abnormal helicase activity 
causes several developmental defects in fungi and plants, such as 
slow growth, premature aging and aberrant assembly of chromo- 
somes [23,24]. These observations suggest a close relationship 
between the abnormalities observed and this helicase. The absence 
of a 1-10 kb sequence from the W4 genome (wild-type) and its 
presence in the telomeric region in the D4 genome (degenerate), 
which displays a high rate of recombination and translocation in 
P. ostreatus, yeast, and humans [25,26,27], supports the hypoth- 
esis that a genomic region harboring a putative helicase was 
transferred from an a strain nucleus to a b strain nucleus to form 
the b^' strain. Thus, this putative helicase might be required for 
vegetative growth and fruiting-body development in F. velutipes. 
Likewise, the transfer of the foreign helicase gene might imbalance 
its activity and result in the observed abnormalities. However we 
cannot rule out the possibility that other region is responsible for 
the degenerations because F. velutipes genome is e (35 
MB)Chromosomal abnormalities can cause the degeneration of 
mushrooms [10,29]. In fungi, extraordinary chromosome rear- 
rangement resulting from insertion/deletion, duplication and 
translocation within and among chromosomes has been reported 
(for review, see [28]). Some of these chromosomal alterations may 
cause fungi to gain a new pathogenicity or adaptation to a novel 
environment [28]. However, in mushroom fungi, in which strain 
stability is critical for stable cultivation, genomic rearrangement 
should be avoided. Recombination predominantiy occurs during 
meiosis, but this strain of edible mushrooms likely degenerated 
during serial passaging (mitosis). In further studies, the mechanism 
by which a genomic region from an a group cell was transferred to 
a b group nucleus to form a degenerated b (b"*) strain should be 
addressed. Moreover, F. velutipes produces abundant arthrospores 
in both monokaryotic and dikaryotic mycelia [12] on media used 
in farms. Thus, the nucleus from a group of mycelia can be 
transferred to a diploid strain (Fvl-5) via plasmogamy. We 
attempted to identify the mobile elements that may cause genomic 
rearrangements [29,30] throughout the flanking sequences of the 
putative hehcase gene, but no such element was detected. Somatic 
recombination has been reported in fungi [29] and may be a 
possible explanation for the mechanism by which a genomic 
region from an "a" group strain, corresponding to degeneration, is 
integrated into a "b" strain genome. The results of the present 
study might be helpful for identifying the potential molecular 
targets responsible for the degenerate phenotype. Further studies 
are required to determine whether helicases are direcdy involved 
in the abnormalities of F. velutipes and to confirm the mechanism 
by which a genomic region is transferred from one nucleus to 
another. 

Supporting Information 

Figure SI Alignment of the sequences corresponding 
to the degenerate-specific region fi-om the degenerate 
(D4), wild-type (Wl) and KACC42780 strains. Gene 
prediction was conducted using FGENESH, and the protein 
function and domain analyses were conducted using the UniProt 
(http://www.uniprot.org/) program and the Protein BLAST 
program of NCBI. The exons of the putative helicase gene are 
shown in bold and the start and stop codons are underlined. * 
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indicates the residues that are identical in the three sequences. ▼ 
and arrows indicate the sequence amphfied using the SCAR 
marker primers that discriminated the degenerates. V and arrows 
indicate the genomic sequences corresponding to the HELICc and 
HepA domains. The telomeric repeating units (TTAGGG) are in 
shaded boxes. 
(PDF) 

Figure S2 Deduced map of the genomic region corre- 
sponding to the degenerate-specific sequence of the 
wild-type (Wl), degenerate (D4) andKAGG42780 strains. 

Unknown but deduced sequences are shown in boxes indicated 
with a broken line . The arrow indicates the direction of 
transcription of the putative hehcase gene. The black boxes in 
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